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Abstract. Fuzzy optimization deals with the problem of determining 'opti- 
mal' solutions of an optimization problem when some of the elements that 
appear in the problem are not precise. In real situations it is usual to have 
information, in systems under consideration, that is not exact. This impre- 
cision can be modeled in a fuzzy environment. Zadeh 20 analyzed systems 
of logic that permit truth values between zero and one instead of the classi- 
cal binary true-false logic. In this framework, satisfying a certain condition 
means to evaluate how close are the elements involved to the complete satis- 
faction. Then, each element in a 'fuzzy set' is coupled with a value in [0, 1] 
that represents the membership level to the set. 

In linear programming some or all the elements that describe a problem may 
be considered fuzzy: objective function, right-hand side vector or constraint 
matrix, the notion of optimality (ordering over the feasible solutions), the 
level of satisfaction of the constraints, etc. Moreover, in integer programs, the 
integrality constraints may be seen as fuzzy constraints. 

Here, we present new complexity results about linear integer programming 
where some of its elements are considered fuzzy using previous results on short 
generating functions for solving multiobjective integer programs. 

1. Introduction 

Integer linear programming is a special case of linear programming in which all 
variables are required to take on integer values only. Its importance is due to the 
amount of real- world problems that can be modeled as integer programs. Some of 
the main well-known applications of integer programming concerns the management 
and efficient use of resources: distribution of goods, production scheduling, machine 
sequencing, capital budgeting, etc. Furthermore, from the mathematical viewpoint 
many combinatoric and geometrical problems in graph theory and logic can seen 
as integer programs. For that reason, and for some others (interesting research 
area...), many textbooks and publications are exclusively devoted to the analysis 
of integer programs, either from a theoretical or practical viewpoint. Some of the 
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text books in this area are those by Schrijver[l6]. Nemhauser-Wolsey [15], Sierksma 
[17j . among many others. 

One of the difficulties when modeling a real-world problem by a mathematical 
programming problem is the possible imprecision of the data or the interpretation of 
the constraints. For instance, suppose in a simple knapsack problem that we want 
to maximize the overall benefits of selecting a subset of items, spending "around" 
$b. We can formulate this problem assuring that we are not going to spend more 
than $6 but in the original problem we were determined to spend around $b and 
maybe increasing the budget a little (such as we consider "around" ) one can obtain 
more benefits. Fuzzy programming deals with this fuzziness which causes difficulties 
in modeling. 

Linear programming problems can be fuzzified in many ways. For instance, 
considering that the coefficients of the objective function or constraints are fuzzy 
numbers or introducing vagueness to the inequalities or to the ordering that de- 
termines the maximization/minimization problem. Moreover, in integer programs, 
the integrality of the variables may be considered fuzzy just fixing how close to an 
integer number one considers that a real number is integer. Some of these pos- 
sibilities for mathematical programming problems are described in [18] (Chapter 
4). 

In general, integer linear problems are NP-hard, then the complexity of the 
fuzzyfication of these problems is as hard as crisp integer programming. 

In this paper we present some complexity results for fuzzy integer problems 
that have never been stated before. We give polynomially results in the sense of 
Lenstra in his well-known result about the polynomially of (crisp) integer problems 
in fixed dimension. Using transformations on fuzzy problems to multiobjective 
integer programs we apply generating functions to prove the polynomial complexity 
of some models of fuzzy integer problems. Short rational generating functions were 
initially used by Barvinok [I] as a tool to develop an algorithm for counting the 
number of integer points inside convex polytopes, based in the previous geometrical 
papers by Brion |5J ,Khovanskii and Puhlikov 12], and Lawrence [13]. The main 
idea is encoding those integral points in a rational generating function in as many 
variables as the dimension of the space where the body lives. Actually, Barvinok 
presented a polynomial-time algorithm when the dimension, n, is fixed, to compute 
those functions. 

The paper is organized as follows. In Section [2] we recall some previous notions 
and results about short generating functions of rational polytopes for multiobjective 
integer programming. We present, in Section [3] complexity results for fuzzy integer 
programs where the inequalities are fuzzified and in Section |4] for integer programs 
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with fuzzy coefficients in the objective functions. Finally, in Section[5]we give some 
conclusions about the results presented through this paper. 

2. Short rational generating functions 

Short rational functions were used by Barvinok [T] as a tool to develop an algo- 
rithm for counting the number of integer points inside convex polytopes, based in 
the previous geometrical paper by Brion [5] . The main idea is encoding those inte- 
gral points in a rational function in as many variables as the dimension of the space 
where the body lives. Let P C R d be a given convex polyhedron, the integral points 
may be expressed in a formal sum f(P, z) — z a with a — (a\, . . . , ad) G PnZ d , 
where z a — z" 1 ■ ■ ■ z°^ d . Barvinok's aimed objective was representing that formal 
sum of monomials in the multivariate polynomial ring Z[zi, . . . ,z n ], as a "short" 
sum of rational functions in the same variables. Actually, Barvinok presented a 
polynomial-time algorithm when the dimension, n, is fixed, to compute those func- 
tions. A clear example is the polytope P = [0, N] C R: the long expression of the 
generating function is f(P, z) = X)ilo z *' anc ^ ^ * s eas y t° see that its representation 
as sum of rational functions is the well known formula 1 ~ z _ z + . 

We recall here some results on short rational functions for rational polytopes, 
that we use in our development. For details the interested reader is referred to 

mm- 

Let P = R™ : Ax < b} be a rational polytope in R™. The integer points 
inside P can be encoded in the following "long" sum of monomials: 

f(P;z)= J2 z ° 

where z a — z" 1 ■ • ■ z" n . Then, to re-encode, in polynomial-time for fixed dimension, 
these integer points in a "short" sum of rational functions in the form 

J(P;z) = 5>^ 

i=i 

where / is a polynomial-size indexing set, and where e G {1,-1} and itj, 6 Z n 
for all i and j (Theorem 5.4 in pQ). 

This encoding tool allows us to present algorithms either for counting integer 
points inside polytopes (see pQ) or for solving single (see [5]) and multi-objective 
integer problems (see [31 H])- The following result states that encoding the set of 
non dominated solutions of a multiobjective integer linear problem can be done in 
polynomial time when the dimension is fixed. It will be useful for our development. 

Theorem 2.1 (|4j). Let A € Z mxn , b G Z m , C = (c x , . . . , c fc ) G Z kxn , and assume 
that the number of variables n is fixed. Suppose P = {x G R" : Ax < b, x > 0} 
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is a rational convex polytope in R". Let MOILPa,c{o) the following multiobjective 
linear integer problem 

(MOILPA,c(b)) max (ci x, . . . , Ck x) =: C x 

(l) s.t.xePnz'l 

whose solutions (called nondominated solutions) are those x* G P fl Z™ such that 
there not exist any other y G P PI Z!f im£/i Cj X* < c,y /or aZZ i = 1, . . . , A; wii/i ai 
Zeasi one strict inequality. 

Then, we can encode, in polynomial time, the entire set of nondominated solu- 
tions for MOILPa,c{o) in a short sum of rational functions. 

One useful result that is used for the proof of the above theorem and for the 
results presented in this paper is the one that states that the computation of the 
short generating function of the intersection of two polytopes is doable in polyno- 
mial time, for fixed dimension, using the generating functions of both polytopes 
(Theorem 3.6 in ■ Basically, it uses the Hadamard product of a pair of power 
series. Given gi(z) = (3 m z m and gi{z) — /J 7 m z m , the Hadamard product 

9 = .9i * 92 is the power series 

g(z) = 2J f] m z m where X] m = /3 m j m . 
roeZ" 



3. Integer programs with fuzzy inequalities 

Let P be a rational polytope, A e Z mx ™, b e Z' Tl and c G II 1 . Consider the 
following integer program where some constraints are fuzzified: 

max c x 

{FIP~ c (b)) s.t. Ax <b 

x g p n z™ 

where < means that the inequalities must be "almost satisfied" . The solutions of 
this problem are then pairs (x, n(xj) where a: is a feasible solution to the prob- 
lem max{cx : x £ P fl Z"} and /i(x) the degree of satisfaction of x to the sys- 
tem of inequalities Ax < b. fi is called the membership function of the fuzzy set 
X = {(x,fj,(x)) : x G X}, where in our case X = P fl Z". The only theoretical 
requirements to /x are: 

(1) fi(X) C M + . 

(2) sup x (i(x) < oo. 

Generally, elements with a zero degree of membership are not listed since it is 
considered that they are far to be a "crisp" (ordinary) element. 



ON THE COMPLEXITY OF SOME FUZZY INTEGER PROGRAMS 



5 



If sup x [i(x) — 1 the fuzzy set is called normal, and of course, any non normal 
fuzzy set can be normalized dividing fi by sup x [J>(x). Then, from now on, we 
consider, w.l.o.g., normal fuzzy sets. 

Then, for FIP^ c (b) we are interested in membership functions that measure 
the satisfaction of the constrains or how far is a solution from the crisp system of 
inequalities, i.e., for each inequality we have a membership function /Xj, i = 1, . . . , m 
such that: 

{0 if ai x > bi + £j 
f(x) if h < a t x <b i + e l 
1 if diX < bi 

where Si are a nonnegative real numbers that measures how we are considering that 
a solution does not satisfied the inequalities at all, and / : W l — > [0, 1]. We assume 
that Ei G Q+, for i = 1, . . . , m, and that /(Z™) C Qn [0,1]. These assumptions 
are not too restrictive since it is usual to consider linear membership functions 
(triangular or trapezoidal fuzzy sets, see [21]) whose coefficients are rational. 

Then, the overall membership function to the system of linear constraints Ax < b 
is given by: 

fi(x) — min fii(x). 



FlPjJb) 



such that 



We are interested in maximizing the linear function cx in 
x G PflZ™ and with the maximum value of the membership function. Since on 
R 2 the componentwise ordering is weak, the above problem is equivalent to find- 
ing the nondominated solutions to a biobjective problem with objective functions 
(cx, fi(x)). Then, we deal with the following equivalent biobjective crisp integer 
problem: 

max (cx,n(x)) 
s.t. xePC\Z\ 

where the fuzzy constraints have been substituted by the membership function fi 
in the objective functions. 

By definition of /x, Problem [2] is equivalent to: 

max (cx,minifXi(x)) 
s.t. x G Pnzi 

and then, equivalent to the following biobjective mixed integer program: 

max (cx,z) 

s.t. z < fii(x), i = l,...,m 

x g P n 1\ 
z e [0,11 
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Let us consider fj,i to be rational linear functions for all i = 1, ...,m. Then, 
because the number of feasible solutions in the x variable is finite ( x G Z" and 
P is a polytope) and for each feasible solution x, z is fixed as min^ fj,j (x) , the 
number of possible nondominated solutions of Problem 0] is finite. Furthermore, z 
is rational because i£P and /ii are rational linear functions. Then, the variable 
z can be transformed to an integer variable y = Mz, with M the least common 
multiple of all the denominators that appear in the inverses of all subdeterminants 
of the matrix defined by the linear functions fJ,i{x), i = 1, . . . ,m. Replacing each 
inequality z < pbiix) by y < Mfii(x), we have that the equivalent problem: 

max (cx,y) 

s.t. y<Mfii(x), i=l,...,m 

(5) 

y G [0, M] n Z 

is a biobjective integer linear program that has the same set of solutions but where 
the solutions in y are related with those in z dividing by M . Furthermore, if c 
is generic for the problem maxjea; : x G P H Z n }, Problem [5] has at most M + 1 
nondominated solutions since for each value of y there is exactly one solution in x. 
Then, this problem may be solved solving M + 1 single objective integer problems, 
one for each of the possible values of y. 

It is worth noting that fuzzy integer programming is NP-hard. Indeed, reduction 
comes from crisp integer programming. It is well-known that finding an optimal 
solution of a general integer program, when the dimension is part of the input, is 
NP-hard (see [TB]). Thus, since we can state that x is a an optimal solution to 
Problem ([2|) if and only if there exists y G [0, M] n Z such that (x, -p-) is an optimal 
solution to the fuzzy integer program FIP^ c (b) the conclusion follows. 

This shows that fuzzy integer programming is as hard as crisp integer program- 
ming. Nevertheless, the situation is even harder because crisp integer programming 
in fixed dimension is polynomial (see [14J ) but fuzzy integer programming is equiv- 
alent to bicriteria linear integer programming (see [1]) which is also NP-hard. 
In spite of that, there is a natural (not easy) way to find all the solutions to 



FIP^ c (b) which is based on solving M + 1 crisp integer problems of the form ^fy, 



fixing y = 0, 1, ... , M. However this approach does not ensure polynomiality even 
in fixed dimension. Here the problem comes from M + 1, the number of problems 
to be solved. This figure might be exponential in the input size and therefore even 
solving each subproblem in polynomial time the overall complexity will be only 
pseudopolynomial . 

The best complexity result that we can state is given by the next theorem. 



ON THE COMPLEXITY OF SOME FUZZY INTEGER PROGRAMS 



Theorem 3.1. We can encode in polynomial time, for fixed dimension, the entire 
set of solutions for FIP^ c (b) 



in a short sum of rational functions. 



Proof. Using Barvinok's algorithm (Theorem 5.4 in [I]), compute the following 
generating function in 2n variables: 

(6) f(x 1} x 2 ):= J2 x^x^ 

where P = {((u,y u ), (v,y v )) G x Z n+1 : u,v G P.cu - cv > 0,y u < 

y v and cu+y u — cv — y v > 1}. P is clearly a rational polytope. For fixed u G Z n , the 
second components, (v,y v ), in the monomial of f(x\,X2) represent 

the solutions dominated by (u,y u ). 

Now, for any function (p, let 7Ti )¥ ,,7r2 iV be the projections of (p(xi,x^) onto the 
x\- and ^-variables, respectively. Thus ^2,7(^2) encodes all dominated feasible 
integral vectors (because the degree vectors of the 2; 1 -variables dominate them, 
by construction), and it can be computed from f(x\,x-2) in polynomial time by 
Theorem 1.7 in [1]. 

Let V(P) be the set of extreme points of the polytope P and choose an inte- 
ger R > max{vi : v G V(P),i — l,...,n} (we can find such an integer R via 
linear programming). For this positive integer, R, and M as described above, let 
r((x,z), (R,M)) be the rational function for the polytope {{u,y u ) G R™ +1 : «, < 
R,i = 1, . . . , n, and y u < M + 1}, its expression is: 

/ 1 7 m+i \ ™ / I r \ 

r((«, ,),(«, M„ = (_ + _) J] (— + j^x) ■ 

Define f{x\,X2) as above, 7r 2) /(a;i) the projection of / onto the second set of 
variables as a function of the xi-variables and F{x\) the short generating function 
of P. They are computed in polynomial time by Theorem 1.7 and Theorem 5.4 in 
PQ respectively. Compute the following difference: 

h(xx) := F(xx) - W2,f(xx). 

This is the sum over all monomials x^^^ where (it, y u ) G P x [0, M + 1] is 
a nondominatcd solution, since we are deleting, from the total sum of feasible 
solutions, the set of dominated ones. 

This construction gives us a short rational function associated with the sum 
over all monomials with degrees being the nondominated solutions for Problem 
As a consequence, we can compute the number of nondominatcd solutions for the 
problem. The complexity of the entire construction being polynomial since we only 
use polynomial time operations among four short rational functions of polytopes 
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(these operations arc the computation of the short rational expressions for f(x±, x%), 
F{x x ), r((x, z), (R, M)) and 7r a ,/(a:i)). □ 



The above result states that the solution of the fuzzy problem can be encoded 
in a short rational generating function in polynomial time for fixed dimension. 
However, to obtain the explicit list of solutions we should expand, as a Laurent 
series, the rational functions that appear in that expression. 

In the following, we present an efficient procedure to obtain the entire set of so- 
lutions for 
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For that, we concentrate on a different concept of complexity 
that has been already used in the literature for slightly different problems. Comput- 
ing maximal independent sets on graphs is known to be #P-hard ([7]), nevertheless 
there exist algorithms for obtaining these sets which ensure that the number of 
operations necessary to obtain two consecutive solutions of the problem is bounded 
by a polynomial in the problem input size. These algorithms are called polynomial 
delay. Formally, an algorithm is said polynomial delay if the delay, which is the 
maximum computation time between two consecutive outputs, is bounded by a 
polynomial in the input size (|10j). 



In our case, a polynomial delay algorithm, in fixed dimension, for solving FIP^ c (b) 
means that once the first solution is computed, either in polynomial time a next 
fuzzy solution is found or the termination of the algorithm is given as an output. 
Next, we present a polynomial delay algorithm, in fixed dimension, for solving 



FIPXM 
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is shown in 



Let L = max{£7, 1" 1 } where U and I are respectively, the largest and smallest 
element that appear in the description of P as a system of inequalities. The pseu- 
docode of a procedure for obtaining the set of solutions of 
Algorithm [TJ 

The following result states the complexity of this algorithm. 

Theorem 3.2. If the dimension is fixed, Algorithm^ is a polynomial- delay method 
to find all the solutions of . 



FlPxJb) 



Proof. Let consider the multiobjective transformation of 



FIPXM 



By definition, P C [0, L] n . Let h(x) denote the short generating function encod- 



ing the nondominatcd solutions of ([5]), rn(x) = JJ^ [— 



1-ar: 



be the short 



generating function of the hypercube TL = J^J[mj, Mi] C M", with m%, Mi <G Q for 

i=l 

i = 1, ... ,n and counting_integer (H) denote the number of operations needed 
to count the integer points encoded in r-H*h (this number is polynomially bounded, 
when the dimension is fixed by Theorem 1.2 in [1]). 
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Algorithm 1: Binary search algorithm for solving MOILP using SGF. 
Initialization: M = [0, L] n C P. 

Step 1: Let Mi, . . . , M2" be the hypercubes obtained dividing M by the 

midpoints of its edges. 

i = 1 

Step 2: 

repeat 

Count the elements encoded in tmAx) * h(x): n^. This is the number of 
nondominated solutions in the hypercube Mi- 
if riMi = then 

if i < 2" then 
I i<-i + l 

else 

| Go to Step 1 with M. the next hypercube to its predecessor hypercube 
end 
else 

if riMi = 1 then 

| Let x* the unique solution in Mi, ND = ND U {x*} and i <— i + 1 
else 

| Go to Step 1 with M = Mi 
end 
end 
until i <= 2" ; 
output: ND 



The algorithm proceeds on a recursive subdivision of hypercubes (subhyper- 
cubes). Starting from the original hypercube [0, L] n , we subdivide it in 2™ sub- 
hypercubes testing whether they contain nondominated solutions (this process is 
done using the short generating function of the corresponding hypercube, r-n, and 
h(x)). 

The subhypercubes that do not contain nondominated solutions are discarded 
from consideration. 

Each subhypercube that contains at least one nondominated solution is subdi- 
vided further until we are led to a family of elements in this subdivision so that 
each one contains either a unique nondominated solution or it does not contain 
any nondominated solution (this search is done using a depth- first search, see |19j). 
Each nondominated solution is added to ND, our current set of solutions. Then, 
all this family of subhypercubes is removed from further consideration. 

The algorithm repeats this scheme with all the elements in the subdivision until 
all of them have been processed: either fathomed or considered because they contain 
a nondominated solution. 
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The entire process is polynomial delay because we are searching in a binary 
tree with at most nlog(L) levels (due to the binary subdivision of the hyper- 
cube [0, L] n ). At each node the algorithm checks in polynomial time, whether 
the branches at this node contain at least one nondominated solution using the test 
given by counting_integer (H) . Therefore, the algorithm only processes those 
branches where it is ensured a final success (i.e. a nondominated solution will be 
found). Hence, from the last solution found, the overall number of operations until 
a next solution is found or finding a certificate of termination is bounded above by 
0(nlog(L)) x count ing_integer (Ti.) . □ 

Remark 3.1. One may think of using a different approach to enumerate the entire 
set of solutions of , 



FlPxJb) 



that may lead to a simpler polynomial delay algorithm. 



The idea would be to use the equivalence between FIP^ c (b) and solving a series 
of M + 1 crisp IP (see the equivalence above). In fix dimension, it is known that 
solving each of these problems is polynomially doable [14] . Then from one solution 
to the next one, generated in this way, the method would need a polynomial number 
of operations and therefore there would be, at most, a polynomial delay between two 
consecutive solutions found. However, this simple method does not guarantee the 
complete enumeration of the set of solution of j 



FlPxJb) 



since each of these IP may 
have multiple optima. The reader may note that to ensure the entire enumeration 
of the set of optimal solutions the method would need to find all the alternative 
optima of each integer problem and this process is equivalent to enumeration of 
integer points in polyhedra which would lead us again to the starting point. 

The following example illustrates a strategy to parametrically fix the membership 
function of the inequalities in 



Example 3.1. Let A e V 

and 



(with rows a%,... ,a m ) and b — (b\, . . . ,b m ) G Z" 



fi(x) = 



for x = (x\ , . • • , x n ) and Pi,qi £ Z + for i = 1 
Then, the hyperplane of M. n+1 given by x 



n+l 



, , TO. 



is the hyperplane that 



contains the i-th face of {x S K™ : Ax < b} (embedded in the hyperplane of x n+ \ 
) and that forms an angle arcsin 



£4 

i=l 



G [0, ?] with the hyperplane 



x n+ i = 0. 

Defining f(x) = min^ fi(x) we have that for each x S R n : 

(1) If Ax < b, then f(x) > 0. And f{x) = if and only if Ax = b. 

(2) If there exists any i £ {1, . . . , m} such that aix > b, then f{x) < 0. 
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Then, defining 

fi(x)=l l + f(x) if-l<f(x)<0 
{ iff{x)<-\ 

clearly fi is a normal membership function modeling the satisfaction of the system 
of inequalities Ax < b. 

In order to minimize as in Problem [H we can assume that /J.(x) is basically 
1 + f(x), with constraints — 1 < f(x) < since in the points inside the polytope 
{x G R n : Ax < b} are catched by those that maximize 1 + f(x) (or equivalently. 
those that maximize f(x) ) and we are not interested in those points that have /J-(x) 
1 1 . so we can obviate the third case in the definition of jj,{x). 



Example 3.2 ([9]). Consider the following problem: 

max 2x\ + 5x2 
s.t. 2xi — X2 < 9 
2x1+8^2 <31 

Xi,X2 G 

We use the membership functions proposed in Example \ 3.1\ with p\ = p2 = 1 and 
qi = 3, (72 = 4 (using the same membership functions as in Figure [7] shows 

the crisp polytope and the area where the membership function is not zero. 
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Figure 1 . The crisp polytope (P) and its maximum deformation 
by the membership function (P) of Example 13.21 



12-2x 1 +x 2 
3 

35 — 2xi —&X2 
4 



Transforming the problem as in ^ we obtain the following biobjective mixed- 
integer problem: 

max (2a;i + 5x2,2) 
s.t. z < 

(8) :; 

xi,x2 e Z + 
z e [0, 1] 

Figure^ shows the feasible region of the above biobjective problem. The bootom of 
that polytope (z = 0) coincides with the embedding of the crisp original polytope 
and in the top (z = 1) appears the maximum deformation of the crisp polytope by 
the fuzzyfication. 
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Figure 2 . Feasible polytope of Problem [5] 



Now, the least common multiple of all the denominators that appear in the in- 
verses of all subdeterminants of the matrix defined by the membership functions is 
M = 12, then, the problem above is equivalent to the biobjective integer problem: 

max f(xi,X2,y):=(2xi+5x2,y) 
s.t. y < 4 (12 - 2xi + x 2 ) 

(9) y < 3 (35 - 2xi - 8x 2 ) 

X\,X2 S 

y e [o, 12] n z 

The nondominated solutions of Problem ([9]) are x\ = 5,x*> = 3,y* = 3 (f* = 
(25,3);, x\ = i,x* 2 = 3,2/* = 9 (f* = (23,9)) and x\ = 3,^ = 3,y* = 12 
(f* = (21,12)^. Then, the set of solutions of the fuzzy problem, with its respective 
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membership values (obtained dividing y by M) are: 

{[(5,3)/0.25],[(4,3)/0.75],[(3,3)/l]} 

4. Fuzzy integer programs with fuzzy objective coefficients 

In this section we deal with linear integer problems where the coefficients of the 
objective functions are fuzzy numbers: 

. „„ . . max cx 

( FIP ^(b)) + DnVn 

s.t. xePnz^i 

where P is a rational polytope and c = (ci, . . . , c„) is a vector of fuzzy numbers. 
First of all, we describe an equivalent way to handle fuzzy numbers that will be 



useful to tackle FIPA,z(b) Let c be a real fuzzy number, for each a G [0, 1] the 
a-cut of c is the set 

c„ = {c e R : (i c - > a}. 

It is clear that c is totally determined by its set of a-cuts for a € [0,1] since the 
membership function that determines c can be identified with this family of sets. 
Actually, the expression for the membership function, given the set of a-cuts is: 

(j,c(x) = sup min{a,xg cl (a:)} 

a£(0,l] 

where Xc a is the characteristic function of the a-cut of c. 

Furthermore, for each a S (0, 1], c a is a closed interval in R: 

c a =[c<?,c%} ae(0,l] 

Although in general, it is necessary to have the complete set of intervals that de- 
scribes the a-cuts, in many well-known families of fuzzy numbers, we can describe 
the fuzzy number just giving a finite subset of a-cuts. For instance, interval fuzzy 
numbers are totally described using one a-cut and triangular or trapezoidal using 
exactly two a-cuts. 

n 

The objective function of our problem, FIPA.c(b) is cx = CiXi that is a fuzzy 



i=i 

number. Let [cf 1 ,cf 2 ] the a-cut for Cj, a S (0,1], i = l,...,rt. Then, using the 
addition and multiplication of a fuzzy number by an ordinary number (see , 
the a-cut for cx is given by: 

n n 
(cx) a = (y ] Cfi Xi, cf 2 Xi] 
i=l i=l 

Now, evaluating a feasible solution x in the objective function means to compute 
a fuzzy number or cquivalently, its a-cuts. To compare two feasible solutions x and 
y, we have to compare fuzzy numbers (that are actually functions), but orderings 



ON THE COMPLEXITY OF SOME FUZZY INTEGER PROGRAMS 



15 



defined over functions are not total, but partial (it may exists functions that are not 
comparable). However, by means of a cuts, we can compare both fuzzy numbers 
(the evaluation by the objective function in x and y) using the equivalent way to 
treat them. 

Let c an d be two real fuzzy numbers, we say that c < d if c" < d" and C2 < d% 
for all a £ (0, 1]. Then, to compare two fuzzy numbers by this partial ordering we 
just need the extreme points of each of its a-cuts intervals, i.e., to compare c and 
d we only need to compare by the componentwise order in R 2 the set of vectors 
{(cf, c% ) : a £ (0, 1]} and {(df , ) : a £ (0, 1]} for each a-cut. 

For our particular case, for two feasible solutions, x and y, we need to compare 

n n n n 

c?i %i and ^ eg y h and ^ c? 2 x, and ^ cf 2 yt. 

i—1 4—1 i—1 i— 1 

Then, we can transform | FIPa , g ( b ) to a continuum family of biobjective problems: 

for each a £ (0, 1] and where c" = (c" 1; . . . , c" n ) and c% = (cfu ■ • ■ , c^) are the 
lower and upper extremes of the a-cuts of c, respectively. 

The set of solutions of these problems are the points x £ P PI HI. such that there 
is no y 6 P fl Z™ with c"a; < cfy and C2 a; < c%y for all a £ (0, 1]. 

In [5] the authors propose a way to reduce this continuum (in a) family of 
problems to a discrete one. Let us consider a%, . . . , £ (0, 1], that we call ranking 
system. Then, we can assume that instead of considering the entire interval (0, 1] 
we consider only a representative set of elements of this interval. In practice, many 
problems are completely determined by a finite subset of a-cuts, so we are not 
loosing information assuming this "discretization" . In those cases, where all the 
a-cuts are needed, we can consider an approximation to the corresponding fuzzy 
numbers with as many elements (finite) in the ranking system as we want. 

The transformation above is exact at the nodes a^ of the representation and 
it has some global errors on [0,1]. Further, it is easy to control the error by 
introducing additional nodes into the representation or by using a sufficiently high 
number of nodes with maxja^ — a^-i} sufficiently small. To control the error of 

i 

the approximation, we can proceed by increasing the number k + 1 of elements in 
the ranking system; a possible strategy is to double the number of points by using 
k = 2 s and by moving automatically to k = 2 S+1 if a better precision is necessary. 

In the standard models, a finite and pre-specified ranking system describes ex- 
actly the fuzzy numbers. For instance, trapezoidal fuzzy numbers (including tri- 
angular fuzzy number as an special case) is totally characterized by two elements 
in the ranking system {ao, 1}. In general it is also true for piecewise linear fuzzy 
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numbers where the ranking system is {ai, . . . , otk — 1}, being a% each one of the 
vertices of the polygonal that gives the membership function. This ranking system 
describes completely the fuzzy number. 

The following example illustrates this idea. 



Example 4.1. Let x be a trapezoidal fuzzy number with membership function given 
by: 



(J&(z) = { 1 








if z < a\ 






z — 




if a± < z 


< 


«2 


0,2- 


-a\ 


1 




if 0-2 < z 


< 


^3 


a i 


— z 


ifa 3 <z 


< 


(7,4 


a.'i - 


-03 







if x > 04 








Figure 3. Trapezoidal membership function. 



Figure shows the membership function of a trapezoidal fuzzy number. Figure [J] 
shows the a^-cut and the 1- cut for some ao S (0,1). From the \-cut, the elements 
ci2 and a% of the fuzzy number are determined. From the otQ-cut, [l\, I2], the equation 
of the line that pass through the points (li, ao) and (02, 1) intersects with the x-axis 
in (oxj 0) and the line that pass through (12, ao) and (03, 1) intersects with the x-axis 
in (04,0). Then, we have completely determined the fuzzy number. 
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Figure 4. ag-cut and 1-cut for a trapezoidal fuzzy number. 



The following result states the complexity of encoding the solutions of FIPA,c(b) 



in a short generating functions when the membership functions involved in the 
problem are piecewise linear functions (or equivalently, that the problem is totally 
described by a finite ranking system) . 

Theorem 4.1. We can encode, in polynomial time, the entire set of solutions 



for FIPa,z{o) in a short sum of rational functions, when a finite ranking system 



describes the fuzzy numbers involved in the problem. 



Proof. Using the a-cuts comparation of feasible solutions, {FIPa,c{o)) is trans- 



formed to the family of problems (FIP^ £ (b)). By hypothesis, it is enough to con 



sider a finite set in the ranking system, so the problem is equivalent to the following 
multiobjective integer problem: 

max v^i ) 2 ' * * * ' 1 ' ^2 ^ ) 

s.t. xePnZ™i 



{FIPIM) 



where the ot\, . . . , is the ranking system. 
Then, the result follows from Theorem 12. II 



□ 



Theorem 4.2. // the dimension is fixed, there exists a polynomial- delay algorithm 



for solving FIPA,z(b) 



Proof. The existence of a polynomial-delay algorithm, similar to Algorithm 1, for 

□ 



FIPA.c{b) follows from the transformation in Theorem 14.2 



In general, fuzzy numbers are not totally described by a finite ranking system 
(this is the case of general LR fuzzy numbers). Then, we propose here a approx- 
imated scheme to solve these problems. Let us consider LR- fuzzy numbers, i.e., 
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fuzzy numbers with the following type of membership function: 

f L(.«i=*.) if z< 01 

where ao < a\ < 02 and L are such that -L(O) = R{0) = 1 and L and i? are strictly 
decreasing continuous functions on [0,1). A example of this type of fuzzy numbers 
is shown in Figure [5] 




Figure 5. LR membership function. 

Introducing as many elements as necessary in the ranking system, we can approx- 
imate the LR fuzzy number above by a continuous piecewise linear fuzzy number. 
Let a* = 4 for some k 6 N \ {0}. Figure [6] shows different choices for the number 
of elements in the system of generators (in the form 041:, with k the number of 
elements). 






Figure 6. Different choices for approximating a LR fuzzy number 
by a polygonal (fc = 3, 7, 15). 



Theorem 4.3. We can encode in polynomial time, for fixed dimension, the entire 
set of solutions for the approximated FIP^ d (b) in a short sum of rational functions. 
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Theorem 4.4. If the dimension is fixed, there exist a polynomial- delay algorithm 



for solving the approximated problem FIP^ £ (b) 



Note that since Theorem 12.11 states that the complexity of multiobjective prob- 
lems does not depends of the number of objective functions (provided finiteness), 
we can increase the number of elements in the ranking system to obtain better 
approximations without increasing the theoretical complexity of the problem. 

The following example illustrates the methodology described above. 

Example 4.2 ([9]). Consider the following problem: 

max ci x\ + 5 X2 
s.t. 2x\ - x 2 < 12 

2x x + 8x 2 < 35 
X\,X2 € Z + 

where C\ is the triangular fuzzy number given by the following membership function: 

f £_! tfl< z <3 

[ otherwise 
Then, the a-cuts for the fuzzy number cx = ci Xi + 5 x% are: 

(cx) a = [(2a + 1) a;i,(5 -2a) 2:1 + 5 aj 2 ] 

that define a the triangular fuzzy number given by the following membership func- 
tion: 

[ otherwise 
Using the ranking system given by {5,1} is enough to solve the problem. After 
transforming our problem to a problem with 4 objective functions, we have that 
Problem (|10[) is equivalent to: 

max (3xi+5x2, 2 x\ + 5 x 2 , 3 x\ + 5 x 2 , 4 x% + 5 X2) 
s.t. 2x\ - x 2 < 12 

2xi + 8x 2 < 35 

Xi, X2 G Z + 

The entire set of nondominated solutions is {(4, 3), (5, 3), (7, 2)}. 



In the following remarks we present extensions of the above problems where 
complexity results can be stated. 
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Remark 4.1. Let us consider the following fuzzy integer problem, where both con- 
strains and objective coefficients are fuzzy numbers: 

max cx 

(12) s.t. Ax <b 

x ePHZf 

Then, 

(1) If the fuzzy numbers involved in ([1211 are totally described by a finite rank- 
ing system, then, the solutions of (|12j) can be encoded in a short generating 
function in polynomial time for fixed dimension. Furthermore, those solu- 
tions can be enumerated using a polynomial delay algorithm. 

(2) // the fuzzy numbers involved in (|12p are not totally described by a finite 
ranking system, then, the solutions of an approximated modification of (| 1 2|) 
(with approximation error as small as desirable) can be encoded in a short 
generating function in polynomial time for fixed dimension. Those solutions 
can be enumerated using a polynomial delay algorithm. 

Proof. The result follows from the following equivalent transformation of Problem 



max (c" x, x, y) 

(13) s.t. xePDZl 

y G [0, M] n Z 

with M as described for Problem ([5]). □ 

Remark 4.2 (Multiobjective Fuzzy Integer Programming). Let us consider the 
following multiobjective fuzzy integer problem, where both constrains and objective 
coefficients are fuzzy numbers: 

max C x 

(14) s.t. Ax <b 

x E P n 7L\ 

where C is a k x m matrix of rational fuzzy numbers. 
Then, 

(1) // the fuzzy numbers involved in (j!4[) are totally described by a finite rank- 
ing system, then, the solutions of (|14[) can be encoded in a short generating 
function in polynomial time for fixed dimension. Furthermore, those solu- 
tions can be enumerated using a polynomial delay algorithm. 

(2) // the fuzzy numbers involved in (fl4|) are not totally described by a finite 
ranking system, then, the solutions of an approximated modification of (| 14[) 
(with approximation error as small as desirable) can be encoded in a short 
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generating function in polynomial time for fixed dimension. Those solutions 
can be enumerated using a polynomial delay algorithm. 

Proof. The result follows from the following equivalent transformation of Problem 

(15) si. iePnZf 

y e [o, M] n z 

where, and c" 2 are the lower and upper extremes of the a-cut of the j-th row 
of C and M as described for Problem ([5]) . □ 

5. Conclusions 

In this paper we present methodologies for solving different models of fuzzy in- 
teger programs analyzing their theoretical complexity. We deal with fuzzy integer 
programs with fuzzy constraints and imprecise costs. The proofs of the results pre- 
sented through this paper are based on the transformations of the fuzzy problems to 
(crisp) multiobjective integer programs and the use of generating functions of ratio- 
nal polytopes. We prove new complexity results about fuzzy integer programming, 
concluding that: (1) Encoding the entire set of optimal solutions of a broad class 
of fuzzy integer programs in a short generating function is doable in polynomial 
time for fixed dimension; and (2) Enumerating these solutions can be done using 
a polynomial-delay algorithm. For problems with imprecise cost where the fuzzy 
numbers involved in the problem arc not totally described by a finite ranking sys- 
tem, we present similar results but considering approximated fuzzy numbers. The 
advantage of the results presented for this approximation is that the theoretical 
complexity does not depend of the global error considered for the approximation, 
and then we can choose approximations as precise as we want. Finally, we give 
similar results also for multiobjective integer fuzzy programs. 
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